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I Abstract. Supersymmetric (SUSY) transformation operators with complex factor- 

ization constants are analyzed as operators acting in the Hilbert space of functions 
i-C . square integrable on the positive semiaxis. Obtained results are applied to Hamil- 

tonians possessing spectral singularities which are non-Hermitian SUSY partners of 
*I3 ' selfadjoint operators. A new regularization procedure for the resolution of the identity 

C$ , operator in terms of continuous biorthonormal set of the non-Hermitian Hamiltonian 

eigenfunctions is proposed. It is also argued that, if the binorm of continuous spectrum 
eigenfunctions is interpreted in the same way as the norm of similar functions in the 
usual Hermitian case, then one can state that the function corresponding to a spectral 
^ ' singularity has zero binorm. 
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Introduction. In many cases the spectrum of non-selfadjoint operators (e.g. 
Hamiltonians) differ from the spectrum of selfadjoint ones by two essential features. 
These are (i) the possible presence of exceptional points and (ii) the possible presence 
of spectral singularities. Note that exceptional points appear already in the finite 
dimensional case (matrices) whereas spectral singularities are a characteristic feature 
of Hamiltonians possessing a continuous spectrum and, hence, they are impossible for 
finite dimensional operators. We think that just by this reason non-selfadjoint operators 
with spectral singularities are studied in much less details. 

Recently, one can notice a growing interest to Hamiltonians possessing a real 
spectrum and spectral singularities [U 121 El H]- Probably it is due to a remark that 
they may produce a resonance-like effect in some experiments [2]. One can notice a 
contradiction in recently published results. In particular, some authors claim that for a 
Hamiltonian possessing a spectral singularity no any resolution of the identity operator 
is possible [HE]- From the other hand there is a thorough analysis of an exactly solvable 
complex potential defined on the real axis where the authors prove that there always 
exists a set of test functions for which the resolution of the identity operator over a 
biorthonormal set of eigenfunctions takes place [3]. 

Note also that, as shown in [5], the method of supersymmetric quantum mechanics 
(SUSY QM, for a review see e.g. |6J) may be very useful for studying both exceptional 
points and spectral singularities since non-Hermitian Hamiltonians possessing these 



SUSY transformations with complex factorization constants 



2 



features may appear as SUSY partners of Hermitian Hamiltonians. Another property 
of non-Hermitian Hamiltonians widely discussed in the current literature is a possible 
real character of their spectrum (for a review see [7]). Below we will consider only non- 
Hermitian Hamiltonians with a real spectrum which are SUSY partners of Hermitian 
Hamiltonians. This possibility is open after the publication of papers [H] where the 
authors showed that SUSY transformations may produce non-Hermitian Hamiltonians 
with real spectra. 

It happens that for a problem on a semiaxis, exceptional points appear if the Jost 
solution with a real momentum is taken as the transformation function (see below). If we 
displace this real momentum from the real axis to the complex plane, the transformation 
function will correspond to a complex factorization constant but the spectral singularity 
disappears from the spectrum of the non-Hermitian Hamiltonian. In such a way one 
may realize a smooth path in the space of parameters, the non-Hermitian Hamiltonian 
H depends on, leading from points where no any spectral singularity is present to a point 
where H possesses a spectral singularity. This opens a way for regularizing the resolution 
of the identity operator expressed in terms of continuous spectrum eigenfunctions of H in 
the point where the spectral singularity is present. To study this possibility in details, 
one has to consider the Hamiltonian in the vicinity of the spectral singularity in the 
space of parameters where no problems with the resolution of the identity operator is 
present and then consider a smooth limit to the singular point. 

To realize the above described procedure one needs to use SUSY transformations 
with complex factorization constants. Therefore we start with a thorough analysis of 
these transformations as operators acting in the Hilbert space of square integrable 
functions defined on the positive semiaxis. 

Below we will show that a singularity inside the integral describing the resolution 
of the identity operator is caused by a vanishing bi-normalization coefficient of 
corresponding continuous spectrum eigenfunction. Therefore, to regularize this 
resolution, it is sufficiently to shift the factorization constant by an infinitesimal value 
from the real axis to the complex plane in the normalization factor of continuous 
spectrum eigenfunctions. Moreover, below we explain what we mean by the bi- 
normalization coefficient of the continuous spectrum eigenfunction and arrive at the 
conclusion that the function corresponding to the spectral singularity has zero binorm. 

SUSY transformations of a real scattering potential with complex factorization 
constants. Let we are given a Hermitian scattering Hamiltonian ho with a real valued 
potential Vq(x) = Vq(x) 



It is initially defined on the domain T>h C C 2 in the space C? := C 2 {1Z + ) of functions 
square integrable on the non- negative semiaxis 1Z + . As the domain T>h one can choose a 
set of finite twice continuously differentiable functions ip(x) with the Dirichlet boundary 
condition at the origin, ip(0) = 0. In this case the closure of h is self-adjoint in C 2 (see, 
e.g., [9]) which we will denote by h = h and its domain by Vt . 




(1) 
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For simplicity we will assume also that h Q has a purely continuous spectrum, 

hoiJ k = E k ^ k , E k = k 2 , k>0. (2) 

Functions i[>k( x ) (generalized eigenfunctions, see e.g. [JO]) satisfy the Dirichlet boundary 
condition at the origin, i[) k (0) = 0, and asymptotic condition at x = oo. We would like 
to emphasize that they are assumed to be real valued, 

r k (x)=M*), k>0. (3) 

Moreover, they form an orthonormal basis (in the sense of distributions) in the space 
C 2 

POO poo 

(ipk\tpk') ■= / dxipl(x)i) k '{x) = (k - k') , / dk\ip k }(ip k \ = 1 (4) 
Jo Jo 

and realize the spectral decomposition of h which is just the colsure h of the differential 
operator h introduced above, 

POD 

h = / dkk 2 \^ k ){^ k \. (5) 
Jo 

Note that ip k (x) £ C 2 but they belong to a wider space of functions defining regular 
functionals over a set of test functions. The action of the operator h is naturally 
extended to the space of these functionals. The interested reader may find more 
information about rigged Hilbert spaces in [TO] . 

Operator ([5]) is defined on a wider domain T)-^ than /i . It consists of all functions 
ip G C 2 such that 

poo 

||M>|| 2 := (M'lM') = / dkk 4 (i)\ip k )(ip k \ip) < oo. 

Operator Hq presents the minimal closed and self-adjoint extension of h . 
Operator (DO) and a non-Hermitian Hamiltonian 

H = -d 2 x + V(x) (6) 

together with its adjoint = —d 2 + V*(x) are related via intertwining relations 

Lh = HL , h tf = L ] H^ . (7) 

Using the real valued character of vq we obtain a complex conjugate form of equations 
© 

L*h = H*L* , h Q L* ] = L* ] H* ] = L* ] H . (8) 

Here we used the property H* = which we assume to hold. We would like to 
emphasize that U here should be understood as formally (Laplace) adjoint to L. Below 
we will describe domains of the operators L, L* and U as differential operators acting 
in the space C 2 and give their closed extensions. 

In the simplest case, that we shall consider below, operator L is a first order 
differential operator 

L = —d x + w (9) 
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with w = w(x) = [\ogu(x)}' being a complex valued superpotential denned with the 
help of a complex valued solution u = u(x) ^ OWx E (0, oo) of the differential equation 
hou = au with a being, in general, a complex factorization constant. The constant a 
participates at the factorization of the Hamiltonians ho and H 

L* ] L = h -a, LL* ] = H-a. (10) 

Complex conjugate form of these relations may also be useful 

L^L* = h -a*, L*L ] = H ] -a*. (11) 

Potential V(x), defining Hamiltonian ([6]), is expressed via the superpotential in the 
usual way, V(x) = vo(x) — 2w'(x). 

Any solution <p> of the differential equation 

H<p E = E<p E , E^a (12) 

may be obtained by the action on a solution of the equation h ip E — Eip E with the 
operator L, cp E — Ltp E . When E = a the function ^ a = 1/m is a solution to equation 
( |T2l) . Hip a = onp a . Another solution to this equation at E = a, <p a , linearly independent 
with ip a , may be obtained from the equation ip a (x)(p' a (x) — (p' a (x)<p a (x) = 1. 

In general, the transformation L violates boundary conditions. The problem is 
simplified essentially if L transforms the eigenfunctions of the Hamiltonian h , ip^, to 
the eigenfunctions of the Hamiltonian H, ipj.. To guaranty this property we will solve 
the boundary value problem defined by the equation ( fl2i) and the boundary condition 

[(f' E {x) + ip E (x)w(x)] x=0 = . (13) 

The initial domain T>h of the operator H as an operator acting in the space C 2 consists 
of a set of finite twice continuously differentiable functions <p(x), satisfying condition 
f|T3|) . The domain V H \ is defined asD^t = (V H )* . Note that W defined on V H \ is not 
adjoint to H in the sense of the inner product. Below we will give their closed extensions 
where this property will take place. 

Equations (TlOT) and ( Till are derived at the level of differential expressions for all the 
operators involved in these expressions. Operators ho, H and have initial domains 
in the space £} described above. Using formulas ( fTUl) and ( TTTT) we can define L, L* and 
L*^ as operators acting in the space C 2 . Evidently, operator L @ is well defined on 
T>h . Therefore we take T>h as an initial domain for L acting in C 2 . Similarly, T>h and 
V H t may be taken as initial domains for L*^ and iJ respectively. 

The functions (p^ = Lip]-, k > solve the equation Hip^ = k 2 ipk, satisfy the 
condition ( fl3|) . and the asymptotic condition at infinity and, hence, they are continuous 
spectrum eigenfunctions for H. Similarly to ipk they do not belong to the space C 2 but 
to a wider space of functions defining regular functionals over a set of test functions. 
Similar to ho the operator H is extended to the space of these functionals. Note also that 
applying L to ipk ^ C 2 we naturally extend the action of the transformation operator 
to the space of functionals. Below we will show that the functions ipk participate at the 
resolution of the identity operator. 
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The function ip Q = 1/u which is a solution to Eq. (TT2"|) . satisfies the condition flTB]) . 
Therefore to have the spectrum of the operator if real for a complex a = —a 2 with 
3(a) 7^ 0, we have to choose the function u{x) such that l/u(x) — > oo as x — > oo. This 
means that a good choice for u(x) is the Jost solution f(k,x) for the Hamiltonian h , 
h f(k,x) = k 2 f(k,x) which has the asymptotics f(k,x) — >■ exp(z/cx) as x — > oo at a 
fixed value of k = —ia, a G C. Thus, 

w(x) = f(—ia,x) — » exp(ax) , a = —a 2 , x — > oo (14) 

with a = d + ib, b and d < 0. For d = 0, it(x) oscillates at infinity, i.e., it 
satisfies the asymptotic boundary condition corresponding to a continuous spectrum 
eigenfunction. As it will be shown below, just in this case a spectral singularity appears 
in the spectrum of if. For d < 0, u(x) — > as x — > oo and the Hamiltonian if has no 
spectral singularities. In both cases the operator H has a purely continuous spectrum, 
E = Ek = k 2 , k > 0. Below we will distinguish the case when H has a spectral 
singularity from that when its spectrum is regular. 

Regular case. Consider first the case d < 0. As it was already mentioned, 
the functions Lip^, satisfy both the equation (j!2p and the boundary condition f|T3|) . 
Therefore for d < the functions 

ip h = N^lAl) h , N k := (k 2 - a) 1 / 2 (15) 

form a continuous biorthonormal basis (in the sense of distributions) in £ 2 , 



\ dx(p k (x)(f k/ (x) = (k - k') , / 
Jo Jo 



((fl>\(fk) = / dxcp k {x)(p k >(x) = (k - k') , / dk\ip k )(ip* k \ = 1 . (16) 



The first property follows from the first factorization relation in ( jTUi) and the eigenvalue 
equation ([2]). The second property is a characteristic feature of a non-Hermitian 
Hamiltonian with a purely continuous spectrum without spectral singularities (see e.g. 
[TT| fl~2]). It may be justified by the following chain of equations 

f°°dk \<p k )(<pi\ = l f°°dk MM (L *)t = L{ho _ a)-\L*y 

Jo Jo K — a 

= (if - a)~ 1 L (f *) j = (ff - a)-\H - a) = 1 . (17) 

Here we used equations ( ITBI) . (j3D, spectral decomposition of the resolvent of ho, first 
equation ([7]) and second equation (TTOl) . 

Applying operator L*' to (fT5l) . using the factorization property (TIP]) once again, 
and eigenvalue equation (ED, one obtains the transformation inverse to (TT5l) 

^ k = N^L^(p k . (18) 

Using the basis ( JTBT) . we obtain the spectral decomposition of if which we will denote 
by if 



oo 

2 



ff = / dfc (19) 



o 
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From here one finds the operator adjoint to H 

POO 

H ] = / dkk 2 \cpl){cp k \. (20) 

Jo 

Operator H is defined on a wider domain than H . It consists of all functions ip G C 2 
such that 

POO POO 

|2 _ /T7, iTj,„\ _ / j/„ l2 / J////2, 



||ify|| 2 = (W^> = / dfcA; 2 / dfc'^(^|^,)(^|^ fc )(^|^)<cx). 

Jo Jo 

Evidently, operator if is densely defined by the relation ( 1T9|) . Therefore if is 
closed. Moreover, since H = H this means that H is closed also (see e.g. [T3]). 
Since = if^fc = k 2 (p k an d -^Vfc = -^V^fc = k 2 (p k , operators H and H are closed 
extensions of differential operators H and if ^ respectively. 

Closure of transformation operators in the space C 2 in the regular case. In the 
space C 2 we have two bases. The functions ip k form an orthogonal basis and tp k form a 
biorthogonal one (in the sense of distributions). Transformation from one basis to the 
other is realized by operators 



U= dk\<p k )(ip k \, Uip k = (p k (21) 
Jo 

and 

POO 

U= / dk\iP k )(cpl\, Uip k = ?p k . (22) 



They have the following property UU = 1, UU = 1. The first relation here is the 
usual resolution of the identity operator over the orthonormal basis ip k (see the second 
equation in (j3])) whereas the second relation is nothing but another form of the second 
equation in ( TT6T) . 

Using these bases one can construct operators 



3 POO 

dkN k \<p k )(i> k \, r= / dkN* k \^ k )(ip k \ (23) 



L* f = / dkN k \i; k )(v* k \, L~ = I dkNZ\<p* k )(if) h \ . (24) 



and 



Moreover, if we introduce the shifted versions of the Hamiltonians h Q and H, 

POO POO 

g -=h -a= dk N k \ip k ) (ip k \ , g 1 :=H-a= dk N k \ip k ) (ip* k \ (25) 
Jo Jo 

and their square roots 

POO POO 

gl /2 = dkN k \ip k ){ip k \, g\ /2 = dkN k \ Vk ){ V * k \. (26) 



then from (I2T1). fl26]). (j23l and from (j22jl. fl25]). pi) it follows that 

y rr 1/2 y*t f r 1/2 

In case when H is Hermitian, these relations reduce to polar decompositions of the 
transformation operators obtained in 
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From (123 p we see that the domain of definition of L consists of all functions ip G C 2 
such that 

roc /»oo 

\\m 2 = dkN k dk'N* k ,(^ k ,)(<p k ,\<p k )(^) <oo. 
Jo Jo 

Similarly, from (J21]) it follows that the domain of definition of L consists of all functions 
ip G C 2 such that 

||L*V|| 2 = / rfA;|iV fc | 2 |(^|^)r<oo. 
Jo 

Since L = L these operators are closed (see e.g. |13j). Furthermore, it is easy to check 
that 

L^ k = Lip k = N k ip k , L*V* = £*V* = N k tp k (27) 

meaning that they are closed extensions of the differential transformation operators L 
an L* f . 

Spectral singularity. The eigenf unctions i[) k , k > of h can be expressed via the 
Jost solution f(k,x) as (see e.g. [15]) 

where F(k) = f(k, 0) is the Jost function for h . Since h is self-adjoint, F(k) does not 
vanish for any real k ^ 0, i.e. F(k) ^ VA; G R \ |15j . From here one finds 

^fc(s) = \/f r^yf [F(~k)Lf(k, x) - F(k)Lf(-k, x)] . (28) 

For any factorization constant a = —a 2 , a = d + ib with 6^0 and d < 0, both 
L/(fc, x) 7^ and Lf(—k, x) ^ VA; > 0. But for d = one has a = —a 2 = b 2 > and 
since w(x) = f(—ia, x) = f(b, x), this implies that Lf(b, x) = 0, but Lf(—b, x) ~ l/u(x) 
and hence I/ipk=\b\ ~ l/u(%) — >■ exp(— z6x) as i -) oo. Therefore the function 
V ? a( a; ) — l/u(x) has the asymptotics v ? a( a; ) — ► exp(— i6x) as a; — >■ oo (i.e., this is the 
Jost solution for if at k = —b) and at the same time, as it was already mentioned, it 
satisfies the boundary condition (fT3"|) . By this reason the point E = k 2 = b 2 belongs 
to the continuous spectrum of H and since the corresponding eigenfunction ip k= \b\{x) 
coincides (up to a nonzero constant) with the Jost solution for H, the point E = k 2 = b 2 , 
according to the definition (see e.g. [121 111 E]) 5 is the spectral singularity for H. By this 
reason and by the definition of the Jost function as the value of the Jost solution at the 
origin, the Jost function for H vanishes at k — —b. Furthermore, the resolvent of H 
grows in C 2 norm in the vicinity of k = — b faster than in the absence of the spectral 
singularity [2] and a singularity appears in the resolution of the identity operator at 

k = m m El- 

Another remarkable property of the function Ltp^x) is the value of its binorm. 
First we would like to stress that the notion of norm (or binorm) has a sense if the 
normalization (or corresponding bi-normalization) integral has a finite value which is 
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not the case for the continuous spectrum eigenf unctions. Nevertheless, the usual physical 
terminology attaches a certain significance to the term "normalization of a continuous 
spectrum eigenf unction" . Namely, according to the first equation (jl]), one says that 
the continuous spectrum eigenfunctions are "normalized to the Dirac delta function". 
Note that the Dirac delta is not a real function but a distribution. Thus, the above 
statement about the normalization of the continuous spectrum eigenfunction ip^ Q should 
be understood in the sense of distributions as a functional (k — ko) which, while acting 
on a test function gives the result $(&o)- Just in a similar way, taking into 

account the first formula (116]) . one can say that the continuous spectrum eigenfunction 
<Pk for a Hamiltonian without spectral singularities is "bi-normalized to the Dirac delta 
function" . Let us consider now what happens if the continuous spectrum eigenfunction 
corresponds to a spectral singularity. 

Taking into account factorization ( TIP]) and eigenvalue equation (j2J), one obtains 

((LVfc)*|£Vv> = (iptlL^ Lilly) = (k 2 - a)(k - k') (29) 

where the use of ( TlOl) has been made. Evidently, this quantity vanishes at k 2 = a = b 2 
and the binorm (129]) . depending on the sign of b, has a first order zero at either k = b 
or k = — b times the delta-function singularity which reduces the right hand side of 
( 129]) to zero since the delta-function belongs to a class of slowly increasing generalized 
functions. We have to note that the right (and, hence, the left) hand side of (129]) should 
be understood in the sense of distributions. This means that acting by this functional 
on a test function of k', one obtains a function of k. If the test function is regular at 
k' 2 = a = b 2 , the limit k 2 — > a = b 2 in ( |29l) gives zero. In other words, the limit 
k 2 — > a = b 2 in ( 129]) is zero functional defined on a set of test functions. In particular, 
as the set of test functions one can choose the space of continuous square integrable 
functions Q(k'), k! > 0. Actually, the set of test functions is larger but here we leave 
open the question of how large it is. We have to note that zero functional is regular. 
This means that it may be presented as an ordinary Lebesgue integral with a function 
J^ik') which is zero almost everywhere times a test function. By this reason the value of 
Tik! = \b\) does not affect the zero value of the functional. In particular, this function 
may even be undefined at k! = \b\ (see an example below). Nevertheless, keeping 
the interpretation described above for the normalization of the continuous spectrum 
eigenfunctions, we say that the continuous spectrum eigenfunction corresponding to the 
spectral singularity has zero binorm. The final comment here is that this result is due 
to the vanishing normalization factor N 2 = k 2 — a at k 2 = b 2 if a = b 2 present in ( 1291 ) 
but not to the value of the bi-normalization integral (fT6]) at k — k' — ±6. 

Similar equation has been obtained for a non-Hermitian Hamiltonian defined in the 
Hilbert space of functions square integrable on the whole real line in [3] (see Eq. (11) of 
that paper). This permits us to hope that the the continuous spectrum eigenfunction 
corresponding to a spectral singularity of a non-Hermitian Hamiltonian has zero binorm 
not only for a class of Hamiltonians considered in this article but also for a more general 
case. 
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Resolution of the identity operator for a Hamiltonian with a spectral singularity. 
For any d < the functions realize the resolution of the identity operator (see the 
second equation in ( fT6l) ) which we rewrite in the coordinate representation 



For d = the integrand contains a first order pole at k 2 = a = b 2 . This is just 
a characteristic feature of the spectral singularity In paper [12] a regularization 
procedure is proposed for guarantying the corresponding resolution of the identity 
operator. In a recent paper [3], using a particular example, a thorough analysis of 
a similar regularization is given. Below, using the just established fact that the origin of 
this pole is the normalization factor for tp^, we propose another regularization for this 
integral. 

For d = (a = b 2 ), instead of the functions ipk given in ( |T5l) . let us consider the 
functions (p>k = (k 2 — b 2 — ze) -1 / 2 !^. with the sign of (infinitesimal) e possibly depending 
on the sign of b (see an example below) which differ from ip^ by the normalization factor 
only and use them in the left hand (and the middle) part of equation f )30|) . The integrand 
in this equation is well defined ^ now, but the right hand side is not the Dirac 
delta anymore. Using the property that this relation should be understood in the sense 
of distributions and the known fact (see e.g. [10]) that there always exists a set of 
test functions such that one can interchange the limit as e — > ±0 with the sign of the 
integral, one restores the necessary behavior of the right hand side of fl30l) taking the 
limit e — » ±0 after the integral is calculated. More precisely, we have first to apply the 
functional at the middle part of (|30|) to a test function and after that take the limit 
e ±0. This statement may be additionally justified by the chain of equations 



obtained in a similar way as the chain f JTTj) . We note that two signs ± before e do not 
mean that both signs are relevant. Only one sign leads to the desired result but it may 
depend on the sign of b. 

We have to note also that the above arguments have no strict mathematical rigor 
and should be considered as an indication that this property may really take place and 
probably may rigorously be proven following the lines similar to that which one uses for 
proving the usual Parseval formula (see e.g. [To]). 

We would like also to mention that the above described procedure is a typical 
method for regularizing divergent integrals (see e.g. [10]) but here we apply it at the 
level of distributions. Similar possibility for regularizing the resolution of the identity 
operator at the singular point is also mentioned in [3] where the authors note that 
the infinitesimal shifting the pole of an integrand from the real axis to the complex 
plane may be considered as a procedure alternative to a deformation of the integration 
contour. 




(30) 




(H-b 2 - ie)- l L(L*) ] = (H-b 2 - ie)~\H - b 2 ) ->■ 1 , e ±0 . 



(31) 
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Exactly solvable example. Let us choose v (x) = 0. The Jost solution for h is 
the simple exponential f(k,x) = exp(ikx) so that u(x) = exp(ax), a G C, a / 0, and 
a = —a 2 . The superpotential is just a constant w = a and the transformed Hamiltonian 
contains the kinetic energy only, H = —d 2 , but the boundary condition at x = for 
equation ( fi~2l) contains a complex number, ^'(0) + a<^(0) = 0. We recognize here an 
example of a non-Hermitian Hamiltonian with a spectral singularity first proposed by 
Schwartz [16] (see also [1]). The functions ifk have the form 



which vanishes at k 2 = —a 2 = b 2 V/c' > if d = 0. Thus, this is zero functional acting 
on a set of test functions of the variable k! . The continuous spectrum eigenfunction 
for k = \b\ is 1/u = exp(-ibx). If we tried to calculate the binorm of 1/u using the 
first equation (fT6|) for k = k' = \b\ we will find that this quantity does not exist since 
the integral of exp(— 2ibx) over the interval (0, oo) does not converge. This is just the 
illustration of the property discussed above that the function defining zero functional 
in the form of a usual Lebesgue integral is undefined for k! = \b\. Nevertheless, we say 
that the continuous spectrum eigenfunction corresponding to the spectral singularity 
has zero binorm in the sense described above. 

Let us calculate now the integral in the left hand side of ( J30l) for a = b 2 > 0. For this 
purpose we shift the parameter b 2 in the denominator of the function (132]) from the real 
axis to either the upper half or the lower half of the complex plane by an infinitesimal 
value e, b 2 — > b 2 +ie. We will see that the sign of e depends on the sign of b. After some 
simple manipulations with trigonometric functions and up to a term proportional to e 2 
in the denominator, one obtains for the integrand in ( 13~TT) the following approximation 



After integrating this expression with respect to k over the interval (0, oo), the first 
term gives just the desired Dirac delta function. Thus, it remains to show that all other 
terms give a vanishing result as either e — > +0 or e — > —0. Not that all integrals here 
are standard provided be > and arguments of trigonometric functions are positive (see 
formulas (3.723.2) and (3.723.3) in [TTj ) . These integrals can easily be continued to the 
case be < and a negative argument of cosine function, 




(32) 




x b 2 cos k{x + y) — ibk sin k{x + y) + —(cos k[x + y) + cos k[x — y)) 



(33) 




c> . (34) 
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dkip k (x)ipk{y) ~ (x - y) + 



±& 2 _ ihz ± !£ 

2 



IF F 
± ^ Tz\x-y\ z . = t_ %h (35) 

2^ 26 v ' 



Here 3a = and the upper sign corresponds to Sf^O and the lower sign to 9ftj0. Thus 
the desired integral may be approximated as 

Z 

with the upper sign for be > and the lower sign for be < 0. Replacing the parameter 
z in the middle term of the right hand side of this equation by its value, one can see 
that this term is equal to zero only for be > 0. For be < it survives in the limit e — > 0. 
The last term is a regular functional acting on a test function <&(?/), y > 0, represented 
by the usual Lebesgue integral. Splitting the integration interval in this integral with 
be > in two parts, (0,x) and (x, oo), one obtains 



%e 
2z 



fX POO 

e -zx / e*»$( y )d y + e w / e - zy <S>(y)dy 

JO Jx 



> 0. 



From here it follows that this is a decreasing function of x provided $(?/)—>• as y — > oo 
and therefore for e — > this term vanishes. Thus, in the limit e — > the right hand 
side of formula (1351) reduces to the Dirac delta for a set of continuous test functions 
vanishing at infinity and the sign of e coinciding with the sign of 6. 

Conclusion. In this paper a careful analysis of SUSY transformations with 
complex factorization constants as operator acting in the Hilbert space of square 
integrable functions defined on the positive semiaxis is given. Obtained results 
are applied to non-Hermitian Hamiltonians which are SUSY partners of self-adjoint 
operators. Using the same arguments, which permit to say that the continuous spectrum 
eigenfunctions of a Hermitian Hamiltonian are normalized to the Dirac delta function, 
we arrived at the conclusion, that the continuous spectrum eigenfunction of a non- 
Hermitian Hamiltonian corresponding to the spectral singularity has zero binorm. 

It is also shown that a singularity inside the integral describing the resolution of 
the identity operator is caused by the vanishing bi-normalization factor (k 2 — a) of 
corresponding continuous spectrum eigenfunction with a = b 2 > at the singular point 
k 2 = b 2 . A new procedure for regularizing the resolution of the identity operator at 
the singular point is proposed. It consists of an infinitesimal shifting of the SUSY 
factorization constant from the real axis to the complex plane in the normalization 
factor of the corresponding scattering function. Obtained results are illustrated by the 
simplest example of a non-Hermitian Hamiltonian having the kinetic energy only but 
a complex boundary condition at the origin. Our results, thus, support the conclusion 
made in [3] that the resolution of the identity operator for a non-Hermitian Hamiltonian 
possessing a spectral singularity may be constructed in contrast to the statement made 
in [HE]. 

Despite this positive result, we have to note that, since, as it has been noticed 
by Lyance [12] . there is no possibility for constructing a Hermitian operator equivalent 
to a non-Hermitian Hamiltonian with a spectral singularity, the relevance of such non- 
Hermitian Hamiltonians as fundamental objects in quantum mechanics is still doubtful. 
Nevertheless, for a non-Hermitian Hamiltonian nearly having a spectral singularity (and 



SUSY transformations with complex factorization constants 



12 



having no exceptional points) such a Hermitian operator exists and can be constructed. 
The scattering cross section for this Hamiltonian has the ordinary physical meaning 
and it may exhibit a resonance behavior near the singular point of its non-Hermitian 
counterpart [18]. 
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